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for Pb + Pb at yfs = 2.76 TeV. We find that for the two different beam energies considered the 
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I. INTRODUCTION 

The experimental results accumulated in these years 
in the ultra relativistic heavy ion collisions (uRHICs) 
first in the experiments conducted at RHIC and more 
recently at the LHC has shown that the elliptic flow 
v 2 = (cos(2 ip p )) = {{pl- P 2 y) / (pl + P 2 y)) , is the largest 
ever observed in HIC [l, 2|. The elliptic flow is a mea¬ 
surement of the momentum anisotropy of the emitted 
particles and it is an observable that encodes information 
about the transport properties of the matter created in 
these collisions. Theoretical calculations within viscous 
hydrodynamics ill and in the recent years also calcu¬ 
lation performed within transport approach have 

shown that this large value of V 2 is consistent with a 
matter with a very low shear viscosity to entropy den¬ 
sity ratio 77 /s close to the conjectured lower bound for a 
strongly interacting system, 77 /s = 1/47T [8j. 

While early studies have been focused on elliptic flow 
generated by the global almond shape of the fireball for 
non central collisions. In the recent years the possibil¬ 
ity to measure experimentally the event-by-event angu¬ 
lar distribution of emitted particle has made possible to 
go beyond such a simplified picture accessing the fluc¬ 
tuating shape that encodes higher order harmonics gen¬ 
erating non zero flows v n = {cos{rup p )) M. Hence 
most of the research activity has been now focused on 
the study of the effects of the fluctuations in the ini¬ 
tial geometry due to the fluctuations of the position of 
the nucleons in the overlap region of the collision [Hh 
E3- Such fluctuations in the initial geometry are sources 
for momentum anisotropies of any n-th order harmonics 
v n = {cos(n<p p )) and in particular of the triangular flow 
n = 3, that especially in ultra-central collisions appears 


as the largest one [lT|, Oil, |]j| . 

The comparison between event-by-event viscous hydro- 
dynamical calculations and the experimental results for 
v n seems to confirm a finite but not too large value of 
Anrj/s ~ 1-3 EEa. However, small values of 77 /s is 
not an evidence of the creation of a QGP phase. A phe¬ 
nomenological estimation of its temperature dependence 
could give information if the matter created in these col¬ 
lisions undergoes a phase transition |2til - [2^ | . Information 
about a temperature dependence of 77 /s can be achieved 
studying the V2{pt) and the high order harmonic v u {pt) 
in a wider range of energies. Similar studies have been 
performed using a transport approach but only for the 
elliptic flow in an approach not incorporating event-by¬ 
event fluctuations [22|, |23| . In this paper we extend this 
analysis to high order harmonics studying the role of the 
77 /s on the build up of v u {pt) using for the first time a 
cascade approach with initial state fluctuations. 

There are several theoretical indications that 77/s 
should have a particular behavior with the temperature 
2(1 12 TL 1234271 ] . As an example in Fig 0] it is shown a col¬ 
lection of theoretical results about the temperature de¬ 
pendence of 77/s. FigJT] shows that in general 77/s should 
have a typical behavior of phase transition with a mini¬ 
mum close to the critical temperature Tc (20l - l23| . On one 
hand at low temperature estimates of 77/s in the chiral 
perturbation theory for a meson gas [24 |25I | , have shown 
that in general 77/s is a decreasing function with the tem¬ 
perature, see down-triangles in FigJT] Similar results for 
77/s have been extrapolated from heavy-ion collisions at 
intermediate energies, see HIC-IE diamonds in FigJT] On 
the other hand at higher temperature T > T c 1QCD cal¬ 
culation have shown that in general 77/s becomes an in¬ 
creasing function with the temperature [26|, |29| , see up- 
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FIG. 1: Different parametrizations for p/s as a function 
of the temperature. The orange area refers to the quasi¬ 
particle model predictions for p/s @j. The three different 
lines indicate different possible T dependencies studied in 
this paper. Symbols are as in the legend. See the text for 
more details. 


equation: 


d^f(x,p) = C[f} + S[/ 0 ] (1) 

where C[f] is the Boltzmann-like collision integral. In 
the result shown in this paper we have considered only 
the 2-0-2 processes and for one component system C[f] 
can be written as, 


C[f] = f ( 

J 2 , 1 ', 2 ' 


(/l'/2' - ff2)\M\ S (p + P2 ~Pl> ~P2>) (2) 


where f 2 v 2 , = / II fc =2,i',2'd 3 Pfc/2i?fc(27r) 3 and M de¬ 
notes the transition amplitude for the elastic processes 
which is directly linked to the differential cross section 
|A4| 2 = 167ts(s — 4 M 2 )dcr/dt with s the Mandelstam 
invariant. Numerically we solve the RBT equation us¬ 
ing the so called test particle method and the collision 
integral is solved by using Monte Carlo methods based 
on the stochastic interpretation of transition amplitude 


triangles and circles in Fig[T| but due to the large error 
bars in the 1QCD results for p/s it is not possible to infer 
a clear temperature dependence in the QGP phase. The 
analysis at different energies of V 2 (j>t) and the extension 
to high order harmonics v n {px ) can give further infor¬ 
mation about the T dependence of p/s. In this paper we 
study and discuss the build-up of anisotropic flows v n in 
ultra-relativistic HIC treating the system as a fluid with 
some p/s(T). This is achieved by mean of a transport 
approach with initial state fluctuations. The paper is or¬ 
ganized as follows. In Section [Til we introduce the trans¬ 
port approach at fixed shear viscosity to entropy density 
p/s. In Section Hill we discuss the initial conditions and 
in particular the implementation of the initial state fluc¬ 
tuations in the transport approach. In Section IIV1 we 
study the time evolution of the anisotropic flows ( v n ) 
and the effect of the p/s(T) on the differential v n (pT). 
Finally in Section [V] we study the correlations between 
the initial asymmetry in coordinate space measured by 
the coefficients e n and the final anisotropy in momentum 
space measured by the anisotropic flows (v n ). In this pa¬ 
per we will show results on v n (pT) for n = 2, 3,4 and 5 
for the two different systems Au + Au at yfs = 200 GeV 
and Pb + Pb at yfs = 2.76 TeV at different centralities. 


II. KINETIC APPROACH AT FIXED SHEAR 
VISCOSITY TO ENTROPY DENSITY RATIO 

In this work we employ the kinetic transport theory to 
study the evolution of the fireball created in relativistic 
heavy-ion collisions. We perform such simulations using 
a relativistic transport code developed in these years to 
perform studies of the dynamics of heavy-ion collisions at 
both RHIC and LHC energies i-A 71 3fil [tI . The evolution 
of the phase-space distribution function f(x,p , t) is given 
by solving the Relativistic Boltzmann Transport (RBT) 


In the standard use of the transport theory one fixes 
the microscopical details of the scattering like matrix cle¬ 
ment or cross sections of the processes to study the effect 
of the microscopical details on the observables. This is 
however not our aim we exploit the cross section a to t as 
a tool to determine the p/s of the system. As shown 
in J35| in the hydrodynamic limit observables like V 2 (pt) 
or spectra don’t depend on the microscopic details en¬ 
coded in |A4| 2 . In agreement with the implicit assump¬ 
tion of hydrodynamics. In such an approach it is pos¬ 
sible to study directly the impact of p/s on observables 
like the anisotropic flows v n {px) which is the main focus 
of this paper. Compared with the viscous hydrodynamic 
calculations a kinetic approach at fixed p/s has manly 
two advantages: first in this approach we start from a 
description in terms of f(x,p) instead of starting from 
T IJ ‘ v (x) and it is possible to include initial non equilib¬ 
rium effects (see [32, [33). Second, this approach is not 
based on an ansatz for the viscous corrections for the 
phase-space distribution function 5f with the limitation 
in the transverse momentum range in order to ensure 
that Sf/f « 1. Also this approach provides a tool 
to study the effect of p/s on the observables in a wider 
range of p/s and in transverse momentum pt . Notice 
also that the kinetic freeze-out can be determined self- 
consistently with an increasing p/s(T) that determines 
a smooth switching-off of the scattering rates. A more 
detailed discussion can be found in previous papers, see 
s m, 0 . The disadvantage of the present approach 
is that hadronization has not yet been included. A more 
general disadvantage is that RBT converge to viscous hy¬ 
drodynamics with the relaxation time typical of a kinetic 
theory. However viscous hydrodynamics with relaxation 
times of kinetic theory have been shown to be in quite 
good agreement with experimental data. 

In order to study the dynamical evolution of the fire¬ 
ball with a certain p/s(T) we determine locally in space 
and time the total cross section a to t needed to have the 
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wanted local viscosity. As shown in |3lj the Chapmann- 
Enskog theory correctly describes the relation between 
rj o T,cr(0),p providing a good agreement with the 
results obtained using the Green-Kubo correlator. In 
the Chapmann-Enskog theory and for a pQCD inspired 
cross sectionWypically used in parton cascade approaches 
SSHUmI, da/dt ~ a 2 /{t — m 2 D ) 2 , the ij/s is given 
by the following expression: 


rj/s 


i5 (p> T " = 


1 (P) 

15 g{a)a to tP ’ 


( 3 ) 


where a = m£>/ 2 T, with run being the screening mass 
regulating the angular dependence of the cross section, 
while 77 (a) is the proper function accounting for the per¬ 
tinent relaxation time r ” 1 = g{a)<JtotP associated to the 
shear transport coefficient and it is given by: 


9(a) 



(y 2 +l)K 3 (2y)-yK 2 (2y) 



,( 4 ) 


where K n -s are the Bessel functions and the function 
h relate the transport cross section to the total cross 
section a tr (s) = a tot h(m 2 D /s) with h{Q = 4£(l+£) [(2£+ 
l)ln(l + 1/C) — 2]. 

In order to study the role of the 77/s ratio and its tem¬ 
perature dependence we consider three different cases: 
one with a constant Airy/s = 1 during all the evolu¬ 
tion of the system dashed line in FigJT] another one with 
Airg/s = 1 at higher temperature in the QGP phase and 
an increasing 77/s in the cross over region towards the 
estimated value for hadronic matter Atttj/s « 6 [ 25 i . [4Q] 
shown by solid line in FigJT] Such an increase of 77/s 
in the cross over region 0.8Tc < T < 1.2 Tq allows for a 
smooth realistic realization of the kinetic freeze-out. This 
is because at lower temperature, according to the formula 
Eq.® a oc (77 /s) _1 i.e. the increase of 77/s towards the 
estimated value for the hadronic matter implies the total 
cross section decrease and this permits to achieve in a 
self-consistent way the kinetic freeze-out. In the follow¬ 
ing discussion the term f.o. means to take into account 
the increase of 77/s at low temperature. The third one is 
shown in FigJT] by the dot dashed line. In this case we 
consider the increase of 77/s at higher temperature with a 
linear temperature dependence and a minimum close to 
the critical temperature with a temperature dependence 
similar to that expected from general considerations as 
shown in Figjl] 


III. INITIAL CONDITIONS 

The main novelty in the present paper is the imple¬ 
mentation of initial state fluctuations in a transport cas¬ 
cade approach. We will consider two systems at differ¬ 
ent centralities: Au + Au collisions at yjS nn = 200 GeV 
produced at RHIC and Pb + Pb collisions at yj snn = 
2.76 TeV at LHC. In particular in this section we dis¬ 
cuss the implementation of the initial state fluctuations 


in the above transport approach. In order to generate 
an event by event initial profile we use the Monte-Carlo 
Glauber model. In this model the Woods-Saxon distri¬ 
bution is used to sample randomly the positions of the 
nucleons in the two colliding nucleus A and B. In this 
way a discrete distribution for these nucleons is gener¬ 
ated. We employ the geometrical method to determine 
if the two nucleons one from the nucleus A and the other 
one from the nucleus B are colliding. Within this method 
two nucleons collide each other if the r elative distance in 
the transverse plane is dr < \J <Jnn /tt where ctnn is the 
nucleon-nucleon cross section. In our calculation we have 
used aNN = 4.2 fm 2 for RHIC and &nn = 7.0 fm 2 for 
LHC. N co u and N part are given by counting the number 
of collisions and the number of participating nucleons for 
each event. The next step is the conversion of the dis¬ 
crete distribution for the nucleons into a smooth one by 
assuming for each nucleon a gaussian distribution cen¬ 
tered in the nucleon position. In our model we choose to 
convert the information of the nucleon distribution into 
the density in the transverse plane pr(x, y) which is given 
by the following sum 


Npart 

Pt(x, y) =C ^2 exp 

i— 1 


(x - Xj) 2 + (y~ Vi) 2 

2c -ly 


( 5 ) 


where C is an overall normalization factor fixed by the 
longitudinal distribution dN/dy while a xy is the Gaus¬ 
sian width which regulates the smearing of the fluctua¬ 
tions and in the following calculations it has been fixed 
to a xy = 0.5 fm. In our calculation we have assumed 
initially a longitudinal boost invariant distribution from 
y = —2.5 to y = 2.5. In the first column of Figj2] it is 
shown the contour plot of the initial transverse density 
at mid rapidity for a given event with impact parame¬ 
ter b = 7.5 fm. The upper panel refers to the system 
Au + Au at yjS nn = 200 GeV and the lower panel to 
Pb + Pb at yj snn = 2.76 TeV. 

The transverse density pr(x,y) fixes the initial 
anisotropy in coordinate space that is quantified in terms 
of the following coefficients e n : 

yj (rj. cos ( nq i )) 2 + (r^ sin ( nq i )) 2 


where xt = \Jx 1 + y 2 and = arctan(y/x) is the polar 
coordinate in the transverse plane. In Fig(5]it is shown 
the initial spatial anisotropies e 2 , £ 3 , £4 and £5 as a func¬ 
tion of the impact parameter. The second coefficient e 2 
shows a stronger dependence with the impact parameter 
with respect to the other coefficients because it acquires 
a contribution due to the global almond shape of the fire¬ 
ball while the other harmonics have most of their origin 
in the fluctuations of the positions of the nucleons. For 
more central collisions b < 2.5 fm the e 2 is even smaller 
than the other harmonics because when the effect of the 
elliptic overlap region disappears it becomes more diffi¬ 
cult to have fluctuations of the positions of the nucleons 
along only one preferential direction. 
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FIG. 2: In the left column it is shown the initial transverse 
density pr{x,y ) at mid rapidity for two typical events for 
Au + Au at y/SNN = 200 GeV (upper panel) and Pb + Pb 
at y/SNN = 2.76 TeV (lower panel). In the right column it 
is shown the corresponding initial temperature in transverse 
plane. These plots are for an impact parameter of b = 7.5 fm. 



FIG. 3: Initial spacial anisotropies e n as a function of the 
impact parameter. Different symbols are for different n. The 
solid lines refer to the Monte Carlo Glauber while the dashed 
ones to the optical Glauber model. 


For the initialization in momentum space at RHIC 
(LHC) energies we have considered for partons with 
transverse momentum pr < po = 2 GeV (3 GeV) a 
thermalized spectrum in the transverse plane. Assum¬ 
ing the local equilibrium the initial local temperature in 
the transverse plane T{x,y) is evaluated by using the 
standard thermodynamical relation pt(x,u) = ^T^jiv 2 


with 7 = 2 x (N 2 - 1) + 2 x 2 x N c x Nf = 40 with 
N c = 3 and Nf = 2. In the right column of Figj2] it 
is shown the corresponding initial local temperature in 
transverse plane. As shown in the central region of the 
fireball for mid peripheral collision we can reach tempera¬ 
ture T « 300 MeV at RHIC and T « 400 MeV at LHC. 
While for partons with pt > po we have assumed the 
spectrum of non-quenched minijets according to stan¬ 
dard NLO-pQCD calculations with a power law shape 
j4li ;42|], In coordinate space the partons with pt > po 
have been distributed according the binary collisions. 
The initial transverse momentum of the particles is dis¬ 
tributed uniformly in the azimuthal angle. We fix the 
initial time of the simulation to tq = 0.6 fm/c for RHIC 
and r 0 = 0.3 fm/c for LHC. 

In the following discussion, we will consider two differ¬ 
ent types of initial conditions. One consisting in a fixed 
initial distribution by using the standard Glauber model 
as used in previous works, see Hamlin. The sec¬ 
ond one consisting of an initial profile changing event by 
event according to the MC Glauber model as discussed 
before. 

In our simulations we have used N event = 500 events 
for each centrality class. This number is enough to get 
solid results for the spectra, differential elliptic flow and 
high order flow coefficients v n {pt). For the study of the 
correlations between the initial e n and the final v n that 
will be shown in the next section we have extended this 
analysis to 10 3 events. The inclusion of the initial state 
fluctuations introduce a further difficulties because in or¬ 
der to get stable results we need to have a good sampling 
of the initial geometry event by event and this is con¬ 
trolled by the total number of test particles N tes t ■ Fur¬ 
thermore an irregular initial profile need a good calcula¬ 
tion grid resolution. We have checked the convergency of 
our results for V 2 , iq and tq with the lattice spacing of the 
calculation grid and N test . We found the convergency for 
a grid with a transverse area of the cell At = 0.12 fm 2 
and Ntest = 2 • 10 6 as total number of test particles per 
event. 

The elliptic flow V 2 (pr) and the high order harmonics 
« 3 (p t ), V 4 (pr) and v$(pr) have been calculated as 

v n = (cos[ra(</>- ’F„)]} (7) 


where the momentum space angles 4G are given by 


4a 


1 

= — arctan 
n 


(sin (rKp)) 
(cos (ntjjj) 


( 8 ) 


In this section first we discuss the comparison between 
the Glauber model w/o initial fluctuations with the MC 
Glauber with fluctuations. Without initial state fluctu¬ 
ations only even harmonics can be generated therefore 
we will consider here only V 2 (pr) and iq {pr)- In the 
left panel of FigQ] we compare the differential elliptic 
flow V 2 (pt) obtained with an initial state that changes 
event by event according to the MC Glauber model (solid 
lines) as discussed in details the previous section with the 
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FIG. 4: Results for Au + Au collisions at snn = 200 GeV 
for mid rapidity. Left: differential elliptic flow V 2 (pr) at mid 
rapidity. The solid lines refer to the case with initial state 
fluctuations the dashed lines are for the case without initial 
fluctuations. Thick lines are for b = 2.5 fm while thin lines 
for b = 7.5 fm. Right: differential va(pt) at mid rapidity 
with the same legend as in the left panel. 


one obtained for the case with an averaged initial profile 
(dashed lines). These results are for Au + Au collisions 
at i/s = 200 GeV and for 20 — 30% centrality class. In 
these calculations we have considered Anri/s = 1 at high 
temperature and an increasing rj/s at lower temperature 
as shown by the red solid line in FigJT] As shown for 
mid peripheral collision (with b = 7.5 fm) the effect of 
the fluctuations in the initial geometry is to reduce the 
V 2 (jpt) of about 15%, despite the same initial eccentric¬ 
ity e 2 in Glauber and MC Glauber, see green solid and 
dashed lines in Fig{3] The reduced efficiency in building 
up the V 2 (pr) is related to the fact that for an irregular 
geometry in the transverse plane the pressure gradients 
generate also a small counter-flow towards the inner part 
of the fireball reducing the azimuthal anisotropy in mo¬ 
mentum space due to the global almond shape. The in¬ 
troduction of the fluctuations in the initial geometry play 
the role to generate the higher order harmonics in partic¬ 
ular the odd harmonics which were absent by symmetry 
in the averaged initial configuration. In the right panel of 
Fig|4]we show the same comparison for the quadrangular 
flow v&( pt)- We observe an opposite behaviour: the ini¬ 
tial state fluctuations increase the final V 4 (pr) by a factor 
of 3. This result is related to the fact that the fluctua¬ 
tions introduce about a factor 3 larger initial £4 as shown 
by the comparison between blue solid and dashed lines in 
Fig El In other words for mid peripheral collisions most of 
V2(pr) comes from the global almond shape while V4(pt) 
comes normally from the initial fluctuations. In fact as 
shown by the black thick solid and dashed lines in the left 
panel of FigElthe effect of the fluctuations is to produce a 
larger V 2 (pt)- From the comparison between thick black 
solid and dashed line in the right panel of Fig El we ob¬ 
serve a non zero V 4 (pr) that was absent in the averaged 



FIG. 5: t n /tn(fo) as a function of the time for Au + Au col¬ 
lisions at yjs nn = 200 GeV (dashed lines) and for Pb + Pb 
collisions at ^/Inn = 2.76 TeV (solid lines). Different sym¬ 
bols refer to different harmonics n. 



FIG. 6: (v„) / {vn)™^ as a function of time at mid rapidity 
and for (20 — 30)% of centrality. Dashed lines are the results 
for Au + Au at y/ SNN = 200 GeV while solid lines Pb -f- Pb at 
y/SNN = 2.76 TeV. Different symbols correspond to different 
harmonics. 

initial profile where the initial £4 ~ 0 (see blue dashed 
line in Fig El). Moreover we observe a low sensitivity of 
v 4 (pt) with the centrality similarly to the experimental 
data at RHIC energies [K)|. Such a behaviour would be 
impossible to explain without initial state fluctuations. 


IV. EFFECTS OF p/s(T) ON THE v n {pr) 

In the first part of this section we discuss the time evo¬ 
lution of the eccentricities e n and the anisotropic flows 
coefficients (v n ) for Au + Au collisions at yjs = 200 GeV 
(solid lines) and for Pb + Pb collisions at 1 /s = 2.76 TeV 
In Fig El we plot the time evolution of the e n normalized 
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FIG. 7: Left panel: time evolution of ( v n ) at mid rapidity re¬ 
spectively for (20 — 30)% centrality collisions. Different sym¬ 
bols are for different harmonics. Right panel: Production 
rate as a function of time at mid rapidity and for the 

same centrality. These results are for Pb + Pb collision at 
\/sjvAr = 2.76 TeV. 


to the initial eccentricity e n (to = 0.6 fm/c) for RHIC 
and e n (to = 0.3 fm/c) for LHC. At very early times the 
small deformation of the fireball in the transverse plane 
decrease linearly with time and at first order of this de¬ 
formation we have that e n oc e n (to) — a n t n ~ 2 . This gives 
the ordering in the time evolution of e n shown in Figj5] 
the time evolution of e n is faster for larger n. 

On contrary v n show an opposite behaviour during the 
early times of the expansion of the fireball. In FigjG] it 
shown the average (v n ) normalized to its maximum value 
at the end of the expansion. The ( v n ) appear later for 
larger n and their development is flatter at early times 
for larger n. Similar results have been obtained in a 
2 + II? transport approach where considerations on the 
early times evolution of the fireball give that (v n ) oc t n 
|43|, |44|. As shown in left panel of FigJ7]we observe that 
in the time evolution of the different harmonics ( v n ) the 
ordering is present also at late times. In right panel of 
Figl7]it is shown the production rate for the different har¬ 
monics and as shown different harmonics have different 
production rates. In particular, we observe that at very 
early time the second harmonic has a non zero value for 
^ 0 at variance with higher harmonics for which 
ss 0. This different behaviour could be the origin 
of the stronger correlation between the final elliptic flow 
Vi and its initial eccentricity £2 that becomes weaker be¬ 
tween the final v n and the initial e n for higher harmonics 
(n > 2), see Section IVl 

Differential flow coefficients v n {pT ) are observables 
that carry out more information about the fireball cre¬ 
ated in the heavy ion collisions in particular because they 
are sensitive to the transport properties of the medium 
like the 77 /s ratio. In the following discussion we will 
study the effect of the 77 /s on the build up of the ellip- 




FIG. 8: Upper panel: differential V 2 (pt) (thick lines) and 
va( pr) (thin lines) green and blue lines respectively at mid 
rapidity and for (20 — 30)% collision centrality. The compar¬ 
ison is between the two systems: Au + Au at yfs = 200 GeV 
(left) and Pb + Pb at yfs = 2.76 TeU(right). The dashed lines 
refer to the case with a constant 77/s = ( 47 r ) _1 during all the 
evolution. The solid lines refer to the case with 77/s = (477)^* 
at higher temperature and with an increasing 77/s ratio at 
lower temperature while the dot dashed lines to the case with 
77/s oc T at higher temperature and with an increasing 77/s 
ratio at lower temperature. Right panel: differential V3(pr) 
in red lines with the same legend as in the upper panel. 


tic flow V 2 (pt) and on the high order harmonics ^(pt), 
v 4 (pt) and v$(pt)- With v u (pt) we mean the root mean 
square \J(v 2 ) as it has been done in experimental data 
using the event plane method. In the upper panel of 
Fig ] 5 ] it is shown the elliptic flow V 2 (pr) (green thick 
lines) and the V 4 (pt) (blue thin lines) at mid rapidity 
and for (20 — 30)% centrality for both RHIC Au + Au 
at yfs = 200 GeV (left panel) and LHC Pb + Pb at 
y/s = 2.7 QTeV (right panel). In general in agreement 
with what has been obtained in viscous hydrodynamical 
calculations, the increase of the viscosity of the medium 
has the effect to reduce both V2 and V4. 

As we can see at RHIC energies comparing the thick 
dashed lines with the solid ones, in the left panel of FigJ51 
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the i >2 (pt) is sensitive to the increase of the 77 /s at lower 
temperature close to the cross over region. In particu¬ 
lar the effect is a reduction of the elliptic flow of about 
17%. A similar trend were observed for the 4-th har¬ 
monic V\ {jpt ) where we have a reduction due to the in¬ 
crease of 77 /s at lower temperature but the effect in this 
case is about a factor two larger then the previous one, 
i.e. about 30 — 40%. The different sensitivity to the 
77 /s can be attributed to their different formation time, 
t v4 > t V2 H3- As shown in Fig IS] each harmonics v n start 
to develop at different times. In particular V 4 has its 
maximum development approximatively at r ~ 3 fm/c 
while the V 2 at r « 1.2 fm/c. This means that differ¬ 
ent harmonics probe mainly different temperatures and 
different value of the 77 /s ratio. Assuming that the first 
few fm/c of the expanding fireball are dominated by the 
ID longitudinal expansion f32| where approximatively 
T(t) = Tc/to/t ) 1 ' 3 we have that when V4 has its maxi¬ 
mum development at about r m 3 fm/c the temperature 
is 1.3T/ at RHIC and 27/ at LHC. In other words this 
tell us that V4 at RHIC energies mainly develops closer 
to the cross over region where 77 /s should increase. 

On the other hand at LHC energies, left panel of Fig| 8 l 
the scenario is different, the elliptic flow is almost unaf¬ 
fected by the increase of rj/s ratio at low temperature (in 
the hadronic phase) as we can see comparing the green 
thick dashed line with the solid one. Instead we observe 
that the increase of 77 /s at lower temperature has a more 
sensitive effect on the V4 (pt) with a reduction of about 
5—10%, see blue solid and dashed lines. Again this differ¬ 
ent sensitivity to the 77 /s in the cross over region between 
V2 and i >4 at LHC are consistent with the results obtained 
at RHIC energies and depends on the different formation 
time of the harmonics in relation to the initial T of the 
system. The greater sensitivity at RHIC energies of both 
V2 and 774 to the 77 /s at low temperature is related to the 
different life time of the fireball. In fact the life time of 
the fireball at LHC is greater than that at RHIC, 8-10 
fm/c at LHC against 4-5 fm/c at RHIC. In general this 
means that at RHIC energies the v n have not enough 
time to fully develop in the QGP phase. While at LHC 
energies we have that the v n develops almost completely 
in the QGP phase and therefore it is less sensitive to the 
dynamics in the cross-over and hadronic region. This 
result were firstly found w/o initial state fluctuation in 
refs.[ 7 i, Hil S 3 but remain similar also with fluctuations. 
The last however allow to study for the first time a sim¬ 
ilar effect also on V3{pt)- 

In Fig[ 8 ] it is shown the effect of an r]/s(T) in the QGP 
phase. In the comparison between the solid lines and 
the dot dashed ones the only difference is in the linear 
temperature dependence of 77 /s oc T for T > Tq while 
at lower temperature we have the same dependence (see 
dot dashed lines in FigO}. As we can see the V 4 at LHC 
is sensitive to the change of 77 /s at higher temperature 
while at RHIC energies the V 4 is completely unaffected 
by this change. In the lower panel of Fig[ 8 ] it is shown 
the triangular flow V3(pt) (red lines) at mid rapidity for 



FIG. 9: Comparison between v n {pr) for mid-peripheral (left 
panel) and central (right panel) collision. Different colors refer 
to different harmonics while solid lines correspond to 47777 / s = 
1 in QGP phase and f.o. and dot dashed lines to 77 /s oc T in 
the QGP phase and f.o. 

(20—30)% centrality and for both RHIC Au+Au at •/s = 
200 GeV (left panel) and LHC Pb+Pb at ^ = 2.76 TeV 
(right panel). In agreement with what has been obtained 
for the even harmonics V2 and V4, we observe at RHIC 
energies a reduction of 773 ( 777 / due to the increase of the 
77 /s at low temperature with a reduction of about 25%, 
while at LHC it is almost insensitive to the change of 77 /s 
in the cross over region. However we observe that at LHC 
the third and fourth harmonics are more sensitive to the 
change of rj/s(T ) with respect to the elliptic flow with a 
deviation of about 10% for V3 and V4 against a less 5% 
for t> 2 . Still it has be noted that such a sensitivity is quite 
small to hoping a determination of the T dependence of 
77 /s from the v u (pt)- 

Very recently it has been possible to access also ex¬ 
perimentally to the ultra-central collisions. The ultra¬ 
central collisions are interesting because the initial e n 
come completely from the fluctuations in the initial ge¬ 
ometry rather than by global geometric overlap region. 
In Fig|9]it is shown the comparison of Vnfpx) produced 
in Pb + Pb at -/s = 2.76 TeV collisions for different 
centralities: left panel for mid-peripheral collisions and 
right panel for central collisions. Different colors are for 
different harmonics. Solid lines refer to the case with 
77 /s = 1/(47 t) in the QGP phase and the increase at low 
temperature as shown in FigJT] by red solid lines while 
the dot-dashed lines refer to the case with 77 /s oc T in 
the QGP phase and the increase at low temperature as 
shown in Fig[l]by blue dot dashed lines. From the com¬ 
parison we observe that at low pr both centralities the 
v n {jpr) are much flatter for larger n. This results is in 
agreement with that obtained in hydrodynamic calcula¬ 
tions where v u (pt) oc p'f [46], On the other hand at 
hight pt for ultra-central collisions we observe that the 








1 



ity while a weaker correlation has been found for higher 
harmonics V 3 ,V 4 with €3 and £ 4 . One explanation for 
the weak correlation observed between V 4 and £4 is that 
for final V 4 there is also a correlation with the initial £ 2 . 
In particular in [47] has been shown that it possible to 
have a good linear correlation between V 4 and a linear 
combination of the initial £2 and £ 4 . 

In this section we discuss these correlations within an 
event by event transport approach with initial state fluc¬ 
tuations. A measure of the linear correlation is given by 
the correlation coefficient C(n , m) given by the following 
expression: 


C{n, to) 


E,(4 - (en))(< - {v m )) 
V^n-^n)) 2 E l «-(v m )) 2 


(9) 


FIG. 10: Ratio between v n (pr) for two different parametriza- 
tions of 77 /s as a function of the order of the harmonic n 
and for pr = 1.5 GeV (open symbols) and pr = 2 GeV (full 
symbols). Solid lines refer to mid peripheral collisions while 
dashed lines to ultra-central collisions. 

elliptic flow V2(pt) shows a saturation while for n > 3 
v ti('Pt) increase linearly with px- This is in qualitative 
agreement with what has been observed experimentally, 
but a quantitative comparison would require the inclu¬ 
sion of hadronization that however would not affect the 
sensitivity to rj/s(T). In particular the sensitivity to the 
value of 77 /s in the QGP phase increase with the increas¬ 
ing the order of the harmonics n in agreement with the 
fact that viscous corrections to v n (px) increase with the 
harmonics [35;]. Furthermore we observe that reduction 
of v n (pT) due to the increase of 77 /s in the QGP phase 
(dot dashed lines) is strongly enhanced for ultra-central 
collisions. As shown in Fia llOl for n > 3 the reduction for 
central collisions is about 30 — 35% against a reduction 
of about 10% for mid peripheral collisions. It is indeed 
remarkable that a 30% effect is determined by a slowly 
linear rising of 77 /s with T as the one considered and de¬ 
picted in Fig[T| In particular in central collisions higher 
harmonics acquire a larger sensitivity to the value of the 
viscosity in the QGP phase. Therefore our study suggest 
that to have information about r//s(T) one should focus 
on ultra-central collisions. This point is further strength¬ 
ened by the study of the correlations between v n and the 
initial eccentricities e n that we discuss in the next section. 


V. CORRELATIONS BETWEEN v n AND e„ 

In recent years, the correlation between integrated V 2 
and high order harmonics V 3 , V 4 with the initial asymme¬ 
try in coordinate space £ 2 , £3 and £4 have been studied 
in the event-b y-ev ent ideal and viscous hydrodynamics 
framework [1 71 l l Of . In general it has been shown that the 
elliptic flow is strongly correlated with initial eccentric- 


where e l n and v l m are the values of e n and v m correspond¬ 
ing to the given event i and evaluated according Eq.s 
© and 0. C(n,m) ~ 1 corresponds to a strong linear 
correlation between the initial e n and the final v m . 

The results shown in this section have been obtained 
with N event = 1000 events for each centrality class a 
total number of test particle per event N tes t = 2 • 10 6 . 
In FigJTTJ it is shown the two-dimensional plots of the 
integrated flow coefficients v n as a function of the cor¬ 
responding initial e„ for each event. The results shown 
are for Au + Au collisions at yj snn = 200 GeV and for 
three different centralities (10 — 20)%, (20 — 30)% and 
(30 — 40)%. The viscosity has been fixed to Anp/s = 1 
plus a kinetic f.o. realized by the increase in r]/s(T) 
as in Figtn As shown in the upper panel we observe 
a stronger linear correlation between £2 and V 2 for mid 
central collisions with a linear correlation coefficient that 
shows a monotonic behaviour with the collision central¬ 
ity from (7(2,2) ss 0.96 for (10 - 20)% to (7(2,2) « 0.89 
for (30 — 40)%. Qualitatively the results are in agree¬ 
ment with the one obtained within a 2+1D viscous hy¬ 
drodynamics, see [49]. In general we observe a slightly 
smaller degree of correlation probably induced by the 
fact that we simulate a 3+1D expansion that can be ex¬ 
pected to contribute to the decorrelation. In the middle 
and lower panel of Fig JIT] we have shown similar plots 
for the third and fourth harmonics. We observe again 
a reduction of the correlation coefficient with the cen¬ 
trality of the collision similarly to V 2 and £ 2 . We obtain 
that the correlation between £3 and V 3 for all the colli¬ 
sion centralities is weaker with respect to that obtained 
for the elliptic flow. Furthermore the fourth harmonic 
flow 1)4 shows a weak correlation with the initial £4 in 
particular for mid-peripheral collisions where the linear 
correlation coefficient is quite weak (7(4,4) < 0.3. Fur¬ 
thermore we observe that the {v n )/e n ratio (see dashed 
lines in Fig fill) decreases when decreases the correlation 
coefficient (7(n, n), i.e. for more peripheral collisions. 

A similar behaviour for the linear correlation coeffi¬ 
cient C(n, n ) is observed at LHC energies for Pb + Pb 
collisions at y/sNN = 2.76 TeV. In Fig [12] it is shown 
the comparison between ultra-central and mid-peripheral 
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FIG. 11: t n and v n for Au + Au collisions at snn = 200 GeV and for three different centrality class. Upper panel: £2 and 
V 2 for (10 — 20)%, (20 — 30)% and (30 — 40)% from left to right. Middle panel: £3 and V3 for the same centralities. Finally in 
the lower panel £4 and 114. In these calculations we have fixed 77 /s = l/( 47 r) for hight temperature and the kinetic f.o. at lower 
temperature (see solid line in Fig[l]). 


collision at LHC energies at tJsnn = 2.76 TeV. In ultra¬ 
central collisions the (v n ) are more correlated to the ini¬ 
tial e n than at peripheral collisions and very interesting 
differences emerge looking also at higher harmonics. 

In order to better visualize and discuss such differ¬ 
ences we have plot in Fig [13] the C(n,n) as a function 
of the impact parameter for both RHIC (dashed lines) 
and LHC energies (solid lines). As shown the linear cor¬ 
relation coefficient is a decreasing function of the impact 
parameter for all the harmonics. However, as we can see 
comparing the dashed and solid lines show, at LHC en¬ 
ergies there is a stronger correlation between e n and v n 
for all n with respect to RHIC energies. We observe that 
V 2 and £2 have the same degree of correlation for both 
RHIC and LHC energies while a lower degree of correla¬ 


tion it is shown for higher harmonics n = 3,4 and n = 5. 
More interesting is the fact that for ultra-central colli¬ 
sions at LHC the linear correlation coefficient C(n , n) re¬ 
mains above 0.9 for n = 2,3,4 and even the n = 5 shows 
large C(n,n) = 0.85. This can be visualized also in the 
right panel of Fig fT^l where the (v n ,e n ) correlation plot is 
shown in mid-peripheral (20 — 30%) collisions (left panel) 
and in ultra-central collisions (0 — 0.2%) (right panels). 

The strong correlation observed for ultra-central col¬ 
lisions means that the value obtained for ( v n ) and its 
dependence with the harmonics n for those collisions is 
strongly related to the value of the initial asymmetry 
measure e n . In particular this could imply that the struc¬ 
ture of the v n (pT ) at LHC where C(n, n) ~ 1 carry out in¬ 
formation about the initial geometry of the fluctuations. 
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FIG. 12: e„ and v n for Pb + Pb collisions at ^/snn = 
2.76 TeV and for (20 — 30)% and (0 — 0.2)% centrality cut 
respectively righ and left panel. In these calculations we have 
fixed r\/s = l/(47r) for hight temperature and the kinetic f.o. 
at lower temperature (see solid line in Fig[T]). 



FIG. 13: Correlation coefficient C(n,n) as a function of the 
impact parameter b. Different symbols refer to different har¬ 
monics n. In particular circles, triangles, squares and dia¬ 
monds refer to n = 2, 3, 4 and 5 respectively. The solid lines 
correspond to Au + Au collisions at y/ snn = 200 GeV while 
dashed lines to the system Pb + Pb at \/snn = 2.76 TeV. 


TABLE I: Linear correlation coefficient C(n, n) for RHIC and 
LHC energies and for different temperature parametritazion 
of 77/s. These results are for (20 — 30)% centrality class. 


C(n, n) 

n 

4nri/s = 1 

Akt]/s = 1 + f.o. 

r//s c xT + f.o. 


2 

0.95 

0.94 

0.93 

RHIC 

3 

0.70 

0.58 

0.65 


4 

0.30 

0.28 

0.31 


2 

0.96 

0.96 

0.96 

LHC 

3 

0.78 

0.78 

0.74 


4 

0.39 

0.38 

0.38 


This joined to the observation that for ultra-central col¬ 
lisions the sensitivity of v n to 77 /s is increased by about a 
factor of 2-3 strongly suggests to focus the experimental 
efforts at LHC highest energy and ultra-central collisions. 

To study the effect of the viscosity and its possible tem¬ 
perature dependence on the correlation we have studied 
how change the correlation coefficient with the different 
parametrizations for 77 /s. In Table Q] we show the re¬ 
sults for C(n, n) for the two energies RHIC and LHC for 
(20 — 30)% centrality class. In general for this centrality 
we observe that at LHC energies and for all the viscosi¬ 
ties considered the degree of correlation between e n and 
v n is greater than the one at RHIC energies. Moreover 
we obtain that at LHC the correlation coefficient is not 
sensitive to the change of the viscosity both at low and 
high temperature. A slight different behaviour we have 
at RHIC energies where the effect of the kinetic freeze 
out is to reduce the degree of correlation between the 
initial e n and the final v n . Furthermore, we have com¬ 
puted the non diagonal components for the linear correla¬ 
tion coefficient C(n,m). We found that (7(2,3) ps 0. and 
(7(3,4) ps 0 for all the range of centralities explored which 
means that there is no linear correlation between V 2 and 
£3 and V 3 and £ 4 . A different behaviour we observe for 
(7(4, 2) which is seen to be an increasing function with the 
centrality (7(4,2) ps 0.02 for central collision (b = 0 /to) 
and about (7(4,2) Pb 0.23 at b = 7.8 fm. This means 
that in more peripheral collisions the 4th harmonic V 4 
has some contamination of £2 and it is not driven only 
by £4 as already suggested in [47|. 

Some interesting properties of the v n distributions can 
be inferred by studying the centrality dependence of the 
relative fluctuations a Vn /(v n ). In FigfT-Hi) it is shown the 
( Npart ) dependence of the ratios cr Vrl /(v n } and cr en /(e„) 
where a Vn and a en are the standard deviation respec¬ 
tively for v n and e n . As shown for n = 2 we observe 
a strong dependence of the relative fluctuations with 
the centrality of the collision with <j v2 /(v 2 ) ~ 0.4 for 
( N P art ) ~ 130. For more central collisions this ratio ap¬ 
proaches the value expecte d for a 2 D Gaussian distribu¬ 
tion where u Vn /( v n) = — 1 ~ 0.523 shown 

by dashed line in Fig[l4] For higher harmonic n = 3,4 
and 5 as shown from Figlldb) to FiglLffl) the values of 
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FIG. 14: From panel a) to d) cr Vn /(v n ) (full symbols) and 
ff e „/(f„) (open symbols) as a function of {N par t) respectively 
for n = 2, 3, 4 and 5. The dashed lines indicate the value 
^/4/7 t — 1 expected for a 2D Gaussian distribution. These 
results are for Pb + Pb collisions at ^ /snn = 2.76 TeV. 


cr Vn /( v n) are approximatively the same of the ones of the 
initial geometry with <J Vn /(v n } « <U n /(e n ) and they are 
almost independent of the collision centrality and for all 
the cen tralities studied they are very close to the the 
value y/^/Tt — 1 shown by the dashed lines. These re¬ 
sults imply that the distributions of 1 ) 3 , 774 and 775 for 
all the centrality range studied are consistent with the 
fluctuation-only scenario discussed in [50] and these fluc¬ 
tuations are related to the fluctuations of the initial ge¬ 
ometry. On the other hand, the distribution of V 2 is 
close to this limit for most central collisions while for 
mid-peripheral collisions there is a contribution coming 
from the global average geometry. 


and at LHC for Pb + Pb collisions at y/s = 2.76 TeV. We 
have found that at RHIC the v u (pt) are more affected 
by the value of 77/s at low temperature (T < 1.2Tc) and 
the sensitivity increases with the order of the harmonics. 
At LHC we get a different effect, all the v n (pT) develop 
in the QGP phase at and are not affected by the value 
of 77/s in the cross-over region. However the sensitiv¬ 
ity to the T dependence of the 77/s is quite weak, more 
specifically a constant 77/s = 0.08 or an 77/s ex T induce 
differences in the V 2 of at most a 5% and of about a 10% 
in 773,774,775. The novel result from our analysis is that 
such a scenario changes for ultra-central collisions where 
found an enhancement of the sensitivity of the v u(pt) 
that for n = 3,4,5 reaches about a 30%. We have also 
studied the correlation between the initial asymmetry in 
coordinate space, measured by e„, and the final asymme¬ 
try in momentum space given by (v n ). We have found 
that larger is the collision energy larger is the degree of 
correlation between e n and (v n ). At LHC there is signifi¬ 
cantly more correlation than at RHIC. For both collision 
energies considered and in all the range of impact pa¬ 
rameter studied the V 2 is strongly correlated with the £2 
with the linear correlation coefficient C(2,2) ss 0.95. The 
degree of correlation between e n and the corresponding 
(t 7„) decrease for higher harmonics. Moreover, in ultra¬ 
central collisions we found that C(ti, n) > 0.9 for n = 2,3 
and 4 which imply that the v n oc e n and they carry out 
the information about the initial geometry of the fireball. 
These results joined with the fact that in ultra central 
collisions the v u (pt ) have a large sensitivity to the 77/s 
ratio strongly suggest to focus the experimental effort to 
these collision centrality where it is possible to get bet¬ 
ter constraint on the value of 77/s in the QGP phase and 
having a new insight on the initial state fluctuations. 
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Using an event-by-event transport approach we have 
investigated the build up of the anisotropic flows v u (pt) 
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